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O : ABSTRACT 

o 

(N 

^JC^ • We investigate whether the magnetic field of an accreting neutron star may be 

diamagnetically screened by the accreted matter. We assume the freshly accumulated 
material is unmagnetized, and calculate the rate at which the intrinsic stellar magnetic 
flux is transported into it via Ohmic diffusion from below. For very high accretion 
, rates M (larger than the Eddington rate ME^d)) Brown and Bildsten have shown that 

00 I the liquid ocean and outer crust of the neutron star are built up on a timescale much 

' shorter than the Ohmic penetration time. We confirm this result and go further to 

P\J ■ calculate the screening, both in this limit and at lower accretion rates, where the 

^ ■ screening declines. Considering the liquid ocean and outer crust, we find that the 

O . Ohmic diffusion and accretion timescales are equal for M f« O.lMEdd- 

^ . We calculate the one-dimensional steady-state magnetic field profiles, and show 

that the magnetic field strength decreases as one moves up through the outer crust 
O I and ocean by n orders of magnitude, where n ~ M/0.02 ME^d- We show that these 

' profiles are unstable to buoyancy instabilities when B ^ lO^'^-lO^^ G in the ocean. 

^ ■ This provides a new limit to the possible strength of any buried field. Most steadily 

accreting neutron stars in low-mass X-ray binaries in our Galaxy accrete at rates where 
^ . screening would be effective if the accretion and magnetic geometry is as we portray. 

However, we cannot definitively state that the magnetic fields of these objects are 
screened at high levels. If screened, then the underlying field will emerge after accretion 
halts, on a timescale of only 100-1000 years, set by the Ohmic diffusion time across 
the outer crust. The definitive statement we can make is that magnetic screening is 
ineffective for M < 10~^MEdd! so that, no matter how the accreted material joins 
onto the star, the underlying stellar field should always be evident. In this respect, we 
point out the only known persistently-pulsing accreting X-ray millisecond pulsar, SAX 
J1808.4— 3658, has an accretion rate of M ~ 10~'^MEdd) f^-r below the realm where 
magnetic screening can play a role. 
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1. 



Introduction 



Accreting neutron stars fall generally into two classes (see White, Nagase, & Parmar 1995). 
Those in high mass X-ray binaries (HMXBs) are X-ray pulsars, with dipole magnetic fields 
B ~ 10^^ G, large enough to disrupt the accretion flow and channel it onto the polar caps. Most 
of the neutron stars in low mass X-ray binaries (LMXBs), however, show no direct evidence for a 
magnetic field {B ^ 10^° G). The accretion disk is believed to extend almost to the stellar surface, 
allowing these neutron stars to be spun up to millisecond periods by accretion. Thus LMXB 
ncTitron stars are favored candidates for the progenitors of millisecond radio pulsars (Bhattacharya 
1995), and as such are expected to have magnetic fields B ~ 10^-10^ G. 

Only one LMXB so far has shown a pulse indicative of magnetically-controlled accretion, the 
newly-discovered transient SAX J1808. 4-3658. This source shows persistent 401 Hz pulsations, 
making it the first accretion-powered millisecond X-ray pulsar (Wijnands &: van der Klis 1998; 
Chakrabarty Sz Morgan 1998). Interpreting the pulsations as disruption of the accretion flow by a 
magnetosphere leads to an estimated magnetic field strength of 10^-10^ G (Psaltis &; Chakrabarty 
1999), suggesting that this source will become a millisecond radio pulsar once accretion ceases 
(Bildsten & Chakrabarty 2001). 

The origin of such low external magnetic fields ^ 10^° G is not understood, but is believed 
to be directly connected with accretion of matter onto the neutron star. This suggestion is 
based on the observations that isolated neutron stars show no evidence for dramatic field decay, 
and that low field radio pulsars occur almost exclusively in binaries (Bhattacharya &: Srinivasan 
1995). Utilizing binary evolution calculations, Taam &; van den Heuvel (1986) studied a number 
of systems, finding that magnetic field strength decreases with the amount of accreted matter, 
although Wijers (1997) argues against this interpretation. A number of different mechanisms 
for achieving this reduction in magnetic field strength have been proposed (see Bhattacharya &; 
Srinivasan 1995 for a review). 

One possibility that could explain the low dipole fields and the lack of pulsations from these 
sources is that the external dipole magnetic field is "buried" or "screened" by the accreted matter. 
Bisnovatyi-Kogan &; Romberg (1974) suggested that pulsars in binaries may have low magnetic 
fields because of this effect. Romani (1990, 1995) outlined a picture in which matter initially 
accretes onto the polar cap, eventually spreads under the weight of the accumulated matter, and 
advects magnetic field inwards. He proposed that the magnetic fiux would be trapped in the crust 
and screened by further accretion of diamagnetic matter. 

In this paper, we test this possibility. We consider accretion of unmagnetized matter onto a 
neutron star with a buried magnetic field, and ask whether the field remains buried, or whether 
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the newly accreted material becomes magnetized faster than the rate at which it is accreted, 
allowing the intrinsic stellar field to "leak out". We study the competition between downwards 
compression and advection by accretion at the local rate, m, and upwards Ohmic diffusion. We 
carry out this calculation for the ocean and outer crust of the neutron star that are in steady-state 
once accretion has been occurring for an extended period. We show that, in steady-state, the 
magnetic field strength decreases through the outer crust and ocean by n m/0.02 mEdd orders of 
magnitude, where rriEdd is the Eddington accretion rate. We show that these profiles are unstable 
to buoyancy instabilities in the ocean when B ^ lO^'^-lO^^ G, giving a new limit to the strength 
of any buried field. 

We find that screening could be effective at the accretion rates of most steadily accreting 
neutron stars. Unfortunately, because of the simplified magnetic field and accretion geometry we 
adopt in this paper, we cannot definitively state that these objects have a screened field. However, 
if so, we show that after accretion ceases, we expect the dipole field to emerge on a timescale of 
100-1000 years, set by the Ohmic diffusion time across the outer crust. It therefore seems unlikely 
that screening alone can explain the low magnetic fields of the millisecond radio pulsars. 

Our robust conclusion is that it is not possible to diamagnetically screen the field in the 
crust for accretion rates ^ 0.02 mEdd- The lack of screening at low accretion rates may explain 
why SAX J1808. 4-3658 is the only LMXB so far to show persistent pulsations. Consideration 
of both the observations and evolutionary state of this binary (Chakrabarty & Morgan 1998; 
Bildsten &; Chakrabarty 2001) makes it clear that the time-averaged accretion rate is unusually 
low, m 10~^ wiEddj below the realm where screening can occur. 

Previous authors have considered the competition between accretion and Ohmic diffusion, but 
in different contexts. Brown &: Bildsten (1998) compared the Ohmic diffusion and advection times 
in the atmosphere, ocean and crust of the polar cap of an accreting X-ray pulsar. In this case, 
the accreted matter arrives on the star already threaded by the magnetic field, and may spread 
laterally if Ohmic diffusion is rapid enough. Considering magnetic field evolution in accreting 
neutron stars, Konar & Bhattacharya (1997) assumed that the currents supporting the magnetic 
field were confined to the crust, and calculated the competing effects of Ohmic dissipation, and 
advection of magnetic flux into denser regions with longer dissipation times (and eventually into 
the superconducting core). 

Another process which acts to transport magnetic flux outwards is thcrmoniagnetic drift (e.g., 
Urpin & Yakovlev 1980; Geppert Sz Urpin 1994). In this paper, we consider the liquid ocean and 
solid crust of the neutron star, for which we flnd that thermomagnetic drift is unimportant; the 
magnetic field profile is set by the competition between Ohmic diffusion and accretion. However, 
in the H/He layer which overlies the ocean, we find that thermomagnetic drift dominates the 
transport of magnetic field. The evolution of the magnetic field in this layer, which periodically 
undergoes a thermonuclear instability giving rise to Type I X-ray bursts, will be discussed in a 
separate paper. 
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We start by describing the thermal structure of the ocean and crust, and what sets the 
electrical conductivity there (§2). In §3, we discuss the evolution of the magnetic field under the 
joint action of accretion and Ohmic diffusion, and in §4 evaluate the characteristic timescales for 
these processes. In §5, we compute steady state magnetic profiles. We show that the magnetic field 
decreases with height in the ocean, by an amount that depends strongly on the local accretion rate, 
rh. In §6, we consider the stability of the steady-state profiles to magnetic buoyancy instabilities. 
We summarize our results and conclude in §7 with the implications for steadily and transiently 
accreting neutron stars. 

2. Microphysics and Structure of the Ocean and Crust 

In this section, we outline the structure of the outer layers of the neutron star, and describe 
how we calculate the electrical conductivity. 

2.1. Thermal Structure 

Most neutron stars in LMXBs accrete hydrogen and helium from their companions at rates 
M 10~^^ — 10~^Mq yr~^. Since we are concerned with the outermost layers of thickness 
^ 100 m, small compared to the radius, we adopt plane parallel coordinates and measure 
the accretion rate as a local rate rh (units g cm~^ s~^). The local Eddington accretion rate 
is rhEdd = 2mpc/(l + X)RaT, where cjt is the Thomson scattering cross-section, nip is the 
proton mass, c is the speed of light, R is the stellar radius, and X is the hydrogen mass 
fraction. We use the Eddington accretion rate for solar composition (X = 0.71) and i? = 10 km, 
rriEdd = 8-8 x 10*^ g cm^^ s^^, as our basic unit for the local accretion rate. For a 10 km neutron 
star, this corresponds to a global rate Medd = 1-7 x 10~^ Mq yr~^. 

The neutron star atmosphere is in hydrostatic balance as the accreted hydrogen and helium 
accumulates, because the downward flow speed {v = rh/p) is much less than the sound speed 
Cg. The pressure varies with height as dP/dz = —pg, where the gravitational acceleration g is 
constant in the thin envelope. We assume in this paper that the magnetic field does not play a 
role in hydrostatic balance, i.e. d{B'^ /87T)/dz <C dP/dz (we check this assumption in §6 when we 
consider buoyancy instabilities). A useful variable is the column depth y (units g cm^^), defined 
by dy = —pdz, giving P = gy. As the matter accumulates, a given fluid element moves to greater 
and greater column depth. We take g^ = g/W^^ cms~^ = 1.9, the Newtonian gravitational 
acceleration for a M = 1.4 M0 and i? = 10 km neutron star. 

The accreted H and He burn at y « 10^-10^ g cm~^ and temperature T Ri 2-5 x 10^ K, 
producing heavy elements which form the underlying ocean and crust. The microphysics of 
these layers has been discussed by Bildsten & Cutler (1995) and Brown &; Bildsten (1998). The 
pressure in the ocean and outer crust is mostly provided by degenerate electrons. For the ions, we 
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incorporate Coulomb corrections to the equation of state using the fit of Farouki &; Hamaguchi 
(1993). The crystallization point is determined by F = (Ze)'^ /kBTa, where a is the interion 
spacing given by 47ra^nj/3 = 1, Ze is the ionic charge, rij is the ion number density (we assume 
only one species of ion is present). For typical conditions at top of the ocean, this is 

where p5 = g cm'^, Tg = T/IO^ K, and the mass of each nucleus is Airip. The ions solidify 

when F = F^ 173 (Farouki &; Hamaguchi 1993), at a depth 



Vc 



2 X 10^^ g cm" 



/F9\ (rrn-^ 



where we assume fuUy-relativistic electrons and use the T = Fermi energy. 

We find the thermal structure of the ocean and crust by integrating the heat equation, 

dT _ 3kF 



(3) 



where F is the outward heat flux, and k is the opacity. The heat flux is set by heat released from 
compression of matter by accretion, and from pycnonuclear reactions and electron captures deep in 
the crust (Brown k, Bildsten 1998). For our models, we assume a constant heat flux F = 150 keV 
per nucleon. The opacity k is calculated as described by Schatz ct al. (1999) and Brown (2000) for 
the crust (see discussion in §2.2). The temperature and pressure at the top of the ocean arc set by 
H/He burning in the upper atmosphere. For accretion rates m < mEdd; the H/He burns unstably, 
and we take the temperature and column depth at the top of the ocean to be the conditions at He 
ignition, as calculated by Gumming &; Bildsten (2000). For rh = mEddj we use the steady-state 
burning model of Schatz et al. (1999), giving y = 10^ g cm~^ and T = 5 x 10^ K at the top of the 
ocean. We take the composition of the ocean and crust to be nuclei with Z = 30 and A = 60. We 
arbitrarily stop our integrations at y = 10^^ g cm~^, thus we include only the outer crust, above 
neutron drip. 

The upper panel of Figure 1 shows the temperature profiles for m = 0.01,0.1 and 1 rfi-Edd- 
For the 0.01 and 0.1 mEdd models, we mark the liquid-solid interface with a dot (the mEdd model 
is still liquid at y = 10^^ g cm~^). In the crust, we show two models for each accretion rate, each 
with a different choice for the thermal conductivity: phonon scattering or electron-ion scattering 
(see §2.2 for discussion of conductivity in the crust). 
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2.2. Electrical Conductivity 

In the relaxation time approximation, the electrical conductivity is cr = neC^ jm^Vc^ or 

. = 7.6 X 10- s- (A) f {^\ 

(e.g., Ziman 1964; Yakovlev & Urpin 1980), where Ue is the electron number density, m^, is the 
effective electron mass, fie is the mean molecular weight per electron, and Vc is the electron 
collision frequency. The lower panel of Figure || shows the conductivity as a function of depth. 

In the liquid ocean, electron-ion collisions set the electrical conductivity. We use the results 
of Yakovlev & Urpin (1980), who give 



_ 8.5 X 10^1 s-i 



where Agj ~ 1 is the Coulomb logarithm, and x = pplrrif-c measures the electron Fermi momentum 
Vf- 

The conductivity in the solid crust depends on the level of impurity. Itoh & Kohyama (1993) 
give the electron-impurity collision frequency as 

VeQ = 1.75 X 10^6 s-1 (I) + ^')'^'' (5) 

where Q = Y1, ^ii^i — {Z))'^/ J2 is the impurity parameter, {Z) = Yl, YiZi/Yl, = ^e/ is the 
average charge (the number abundance of species i \sYi = Xi/Ai, where Xi is the mass fraction), 
and the Coulomb logarithm is order unity. The other contribution to the conductivity is scattering 
with phonons, for which we adopt the results of Baiko & Yakovlev (1995). For temperatures much 
greater than the Debye temperature (T » Qn = 2.4 x 10^ K x^/'^{2Z/A)), they find 



13Aep = 1.24 x 10^8 s-i TgAep, (6) 



where a = 1/137 is the fine structure constant, and Agp is of order unity. The ratio of the impurity 
and phonon scattering frequencies is 

where we choose typical conditions in the crust. 

Previous studies have assumed a low level of impurity, Q <C 1, as might be expected for a 
primordial crust. In this case phonon scattering dominates the conductivity. However, recent 
calculations of the products of H/He burning in the atmosphere show that the impurity level is 
probably much greater. Schatz et al. (1999) calculated the ashes of steady-state H/He burning at 
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local accretion rates rh ^ mEdd, and found (Z) = 24 and Q 100. At such high impurity factors, 
it is not clear whether equation (^) is correct. Brown (2000) suggests that when Q ~ {Z)'^, the 
appropriate conductivity is that for electron-ion scattering in the liquid. 

In this paper, we show results for two limiting cases for the crust: a pure crust with only 
phonon scattering, and an impure crust with electron-ion scattering calculated as for the liquid. 
The conductivity in each case is shown in the lower panel of Figure 1. For each accretion rate, 
we plot the two solutions for the crust: for electron-ion scattering, the conductivity is continuous 
across the liquid-solid boundary; for phonon scattering, the conductivity changes discontinuously.0 
The difference in conductivity between these two cases is roughly a factor of 5. In §4, we show 
that the impurity level has an important effect on magnetic field evolution in the outer crust. 



3. Evolution of the Magnetic Field 

We consider a simplified geometry for the magnetic field, a plane parallel model in which the 
magnetic field B depends only on depth. In this case, the vertical component of the magnetic 
field must be constant since V- S = 0, and we set it equal to zero since we assume the accreted 
matter is unmagnetized. This leaves us with a purely horizontal field, B(z) = B{z)ex (we align 
our X-axis with the horizontal component of the field). 

The magnetic field evolves according to Faraday's law 

^ = -cVxE, (8) 

where the electric field E is given by Ohm's law, 

^ V X B J 

E = + -. 9 

c a 

The electric current is 

c c dB 

J = —VxB = ey — — , 10 

for a given profile B(z), and the downwards flow velocity is given by continuity as 

V = —v{z)ez = —{rn/p)ez (we define v so that it is a positive quantity for the downwards 

flow). Thus E = EyCy, where 

^ vB c dB 
^ c Ana oz 



and equation (^) reduces to its x component 

^ r r>¥^^ 

(12) 



dB _ d 
dt dz 



dB 

vB + r]—- 
oz 



^Recent conductivity calculations of Potekhin et al. (1999) show that the conductivity is continuous across the 
melting point even in this case; however, the older results we adopt in this paper are sufficient for our purposes. 
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where rj = c^/Ana is the magnetic diffusivity. 



It is instructive to rewrite equation (12) as 



D fB\ d ( dB\ , , 

- =7jrUijrh (13) 



Dt \ p J dz \ dz 

where D/Dt = d/dt + t^-V is the advective derivative and we have used the continuity equation 
Dp/Dt = —p'V-v. Equation ( [l^ ) clearly illustrates the different processes acting to change the 
magnitude of B. The left hand side describes advection and compression by the downwards flow; 
the right hand side describes Ohmic diffusion. Without diffusion, B in a given fluid element would 
grow oc p as it was compressed. 



In this paper, we calculate only steady state solutions of equation (|12[). Since we assume that 
the incoming material is unmagnetized, we set Ey = 0, and equation (|l2|) reduces to 

rj^ = -vB. (14) 
dz 

We rewrite this so that pressure is the independent variable; this makes the competition between 
advection and diffusion more transparent. Using hydrostatic balance, and introducing the pressure 
scale height H = y/ p = —{dlnP/dz)~^, equation (|lj) becomes 

dlnB ( H'^\ { m\ t^m 

-- 7 > (15) 



dlnP \ f] / \ y 

where we have defined 



V 

the characteristic Ohmic diffusion time across a scale height, and 



jj2 

tdis = — , (16) 



taccr — . ) (1''') 

m 



the characteristic accretion flow time across a pressure scale height. Equation (|15D shows that, in 
steady-state, tdiff/^accr measures the ratio of pressure scale height to magnetic scale height. For 
example, if tdiflf > ^accD a small magnetic scale height is required for Ohmic diffusion to balance 
accretion. 

The magnetic profile is given by integrating equation (^), since all the quantities on the right 
hand side are known functions of pressure (or height). We find 

B{P) = B{Po)S, (18) 

where Pq is a reference pressure, and 



S = exp 



P 



) 



(19) 
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is the amount by which the magnetic field changes due to screening from the overlying matter. 
We shall refer to the quantity S as the magnetic screening factor. 

The Ohmic diffusion time is independent of B (see Goldreich & Reisenegger 1992 for an 



interesting discussion). Thus equation (15) is linear in B, giving the profile of the magnetic field 
independent of its overall amplitude. In the field burial picture, the amplitude of B in the crust is 
fixed by conserving the magnetic flux in the original dipole field. We do not calculate the details 
of the burial process in this paper; thus we treat the overall normalization of the magnetic field in 
the crust as a free parameter. 

For our assumed field geometry, the Hall term in Ohm's law exactly vanishes. More 
generally, the magnitude of the Hall term depends on the ratio of the electron cyclotron frequency, 
uje = Be/m^,c, to the electron collision frequency, = t^^. For a typical value of Vc, this is 



-er=i ^-^ ^^^ ^ , (20) 




where m* = Ee/c^ is the effective mass of the electron. Typically, uJcT > 1 requires B ^ 10^^ G in 
the ocean and crust; we show in §^ that such large fields are buoyantly unstable. In addition, we 
neglect terms in Ohm's law due to thermomagnetic effects (e.g., Urpin & Yakovlev 1980); these 
are unimportant in the ocean and outer crust. 



4. Accretion vs. Ohmic Diffusion 

In we showed that the steady-state magnetic profile depends directly on the ratio of 
accretion and Ohmic diffusion timescales. In this section, we evaluate this ratio, before presenting 
the steady-state profiles in §^ Figure ^ shows the accretion and diffusion timescales as a function 
of depth, and Figure |3| shows idiff/^accr- We find that the diffusion time is longer than accretion 
time for rh ^ 0.1 rh^dd- We now use simple analytic estimates to understand these results, 
discussing the ocean (j |4.1j ) and crust (§ f4.2[) in turn. 



4.1. Ocean 

First, we evaluate tdis = ^tictH'^ / c? . The conductivity is given by equation (^). The 
degenerate electron pressure is = (meC^)^/(2;)/247r^(/ic)'^ (Clayton 1983), or, in terms of column 
depth y = P/g, 

ys = 3.2 fix) (—) (21) 

where 

fix) = x(2x2 - 3)(1 + x^fl'^ + 3sinh-i x (22) 
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(Clayton 1983), and ys = y/W^ g cm ^. We write the density as p = ^e^emp, giving = 9.7 //g x^- 
The scale height H = y/p is thus 

fix 



H = 170 cm 



giving 

^diff 

The accretion timescale is 
giving 



10^ s 



fix) 



A: 



^diff 



Sgix)A-^ 



1.9 

514 

30 

m 



Lgy 

514/ 

-1 



0.1 rhEdd 



(23) 

(24) 
(25) 



m 



0.1 rh-Edd 
30 



Z 



1.9 

914 



(26) 



where gix) = /(x)/x^(l + x^). Thus, tdis ~ ^accr for m 0.1 rn-Edd, in agreement with Figure |3|. 



In the non-relativistic and relativistic limits, the functions fix) and gix) take the limiting 
values 

8x^/5 X < 1 
2x^ X > 1 



and 



Thus 



fix) 

gix) -- 



8x^/5 X < 1 
2/x X » 1 



(27) 
(28) 



/ ^diff 

V taccT J NR 



A- 



m 



0.1 rh-Edd 



for X <C 1, and 



^diff 

^accr / R 



\ z 



m 



1.9\3/5 



914/ 



(29) 



0.1 mEdd 
30 

'z 



L9y/4 

914/ 



(30) 



for X ^ 1. These analytic estimates roughly agree with our numerical results, and explain the 
scaling of tdiff/^accr with depth; for small x, idiff/^accr increases with depth, whereas for large x, 
tdiff/^accr decreases with depth. The total screening factor S through the ocean is thus insensitive 
to the upper and lower boundaries. The maximum value of tdiff/^accr occurs for intermediate 
values of y, for which the electrons are partially relativistic (x ~ 1, /9 ~ 10^-10^ g cm~'^). 
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4.2. Crust 

We now turn to the solid crust. The behavior of tdiff/^accr is quite different in the crust, 
depending on whether we choose electron-ion scattering or phonon scattering. For electron-ion 
scattering, tdiff/^accr oc g{x) oc 1/x, decreasing with depth, as shown in Figure ^. For phonon 
scattering, 

,,„.2.io=s|!(i|&)(M)VAy, (31) 

Tg \ Aep J \giiJ KfleJ 

giving 

m \ 
0.1 mEdd/ 





3 






1.9\ 


i- 


514/ 





13/, 



ep 



-^ep 



(32) 



Thus tdifr/^accr is constant with depth in this case, in agreement with Figure ^. We see that, as 
suggested by Schatz et al. (1999), tjifr is much smaller in an impure crust (electron-ion scattering) 
than a primordial crust (phonon scattering). 

We do not calculate the conductivity of the inner crust, below neutron drip at y ~ 10^^- 
IqW g cj]2-2_ js^^ these depths. Brown & Bildsten (1998) found that t^is < taccr for m < 5mEdd 
(see their Figure 9 for a continuation of our Figure 3 into the inner crust). Thus the inner crust 
does not give a significant contribution to the screening factor. 



5. Steady-State Profiles 

In section ^, we compared the timescale for Ohmic diffusion with that for accretion. We now 



use the results of §3 to compute the steady-state magnetic profiles, by integrating equation (15). 

Figure ^ shows the result for a range of accretion rates. The solid lines are for models with 
electron-ion scattering in the crust; for 0.01 and O.lmEdd; the dotted lines show models with 
phonon scattering. Rather than fix the overall magnitude of the magnetic field, we plot the ratio of 
the magnetic field at each depth to the magnetic field at the base {y = 10^^ g cm~^). As expected 



from inspection of equation (15), we find that the magnetic and pressure scale heights are equal 



for m ~ 0.1 riiEddi the accretion rate at which tdiff ~ t 



accr 



An analytic estimate of the screening factor for the ocean follows straightforwardly 
from the results of The ratio tdiff/^accr is given by equation (|^). We then substitute 
dlnP = {din f /dx) dx in equation (p!9|), where /(x) is given by equation (p^) and 
df/dx = 8xV(l + x2)V2. As we discussed in ^4.1] , most of the contribution to S comes 
from a; ~ 1. Integrating, we find 
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(l+xf)V2 (l + x2)l/2 



(33) 



where xt (xb) is the value of x at the top (bottom), and we take /ie = 2 and A^i = 1. Note that 
equation (^) gives In S; the factor by which the magnetic field changes through the atmosphere is 
the exponential of this value. 

Using the scalings of equation (|3^), and adopting the prefactor from our numerical results, we 
find that the magnetic field decreases by n orders of magnitude through the ocean, where 

Z 



m 



"""^ 0.02 rhEdd V3oj Vl.gj ■ ^^^^ 
The major contribution to n comes from the ocean. Thus the uncertainty in the conductivity of 
the crust (electron-ion or phonon scattering) changes n by only a small amount, as shown by 
comparing the solid and dotted lines in Figure |^. 



6. Buoyancy Instability 

We now investigate the stability of the steady-state magnetic profiles to buoyancy instabilities. 
We first consider interchange and Parker instabilities in §6T, before including the effects of thermal 
diffusion in 



6.1. Interchange and Parker Instabilities 

The simplest case to consider is the interchange instability in which a magnetic field line 
and associated fluid is lifted vertically, maintaining pressure balance with its surroundings. If the 
new density is less than that of the surrounding fluid, it is buoyantly unstable. Newcomb (1961) 
considered the stability of a horizontal magnetic field in a stratified atmosphere, pointing out that 
it is also important to consider Parker-type modes, in which the magnetic field lifts up, but the 
fluid flows back down the field lines. In the limit of long wavelength along the field (minimizing 
the energy used to bend the magnetic field lines), Newcomb showed that the Parker- type modes 
may be unstable even when the interchange modes are not, although with smaller growth rates. 

Newcomb (1961) showed that the criterion for instability is 

where Pg is the gas pressure. Using the equation of hydrostatic balance, dP/dz = 
d{Pg + B'^/S-K)/dz = —pg, to substitute for —pg, we find instability if (see also Acheson 
1979) 

2 \ / \ f dlnB' 




SirPgl V dz 



-) > N^g, (36) 
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where 

(37) 



N"^ _ d\np 1 d\nPg 



g dz Ti dz 

defines the Brunt Vaisala frequency A^^ (e.g., Hansen & Kawaler 1994). Both thermal buoyancy 
and composition gradients may contribute to A^^ (e.g., Bildsten & Gumming 1998). We assume a 
uniform composition in the ocean, so that A^^ has only a thermal piece, 



,2 9 XT dlnT 



H Xp \ dlny J 

where XQ = din P/dlnQ with the other independent thermodynamic variables held constant. 



We see from equation ( pq) that the instability criterion depends on both how quickly the 
magnetic field decreases with height, and the ratio of magnetic to gas pressure, B"^ /SirPg. To 
understand this, first consider a magnetic field whose strength changes discontinuously with 
height. If the magnetic field above the discontinuity is smaller than that below, pressure continuity 
implies heavy fluid overlying light; a Rayleigh Taylor unstable density jump no matter how small 
B^/SirPg. On the other hand, if magnetic pressure dominates gas pressure, even a small magnetic 
field gradient may create an unstable density profile. 

For the steady-state magnetic profile, the gradient dlnB/dz is given by 
din B/ dlny = tdiff/^accr (eq. IH]). Inserting this gradient into equation (|36|) , we find the 
critical magnetic field strength Be required for instability of the steady-state profile is 

We plot Be as a function of depth in Figure ^ (solid line). We assume P ~ in making this 
plot. The dashed line in Figure ^ shows the steady-state magnetic profile, arbitrarily scaled to 
B = W^^ G at the base. For clarity, we show only models with electron-ion scattering in the crust. 

In the degenerate ocean, Bildsten & Gumming (1998) estimated the thermal buoyancy as 

N^H^ « (40) 



allowing us to make an analytic estimate of Be- We use equation ( |26D for tdifr/^accr, which gives 

(41) 



Br = 7 X 10^° G rJ/^ 



9{x) 



/ ■ X -1/2 



0.01 rn-Edd 



where we take Agj = 1, 514 = 1.9, Fi = 4/3 and ^Ug = 2. For x <C 1, we find B^. is weakly dependent 
on depth. Be oc x^/^ oc y^^^^^ . For 1, Be (x x^, giving 

, \ 1/2 / ■ \ -1/2 

-ill n „ ,1/2 Ms \ W m \ ' 



B, = 4xl0iiGy^r hf [tTT^- (^2) 

* 3 / VO.Ol mEdd/ 



in good agreement with our numerical results. 
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6.2. Doubly-Diffusive Instability 



On small enough wavelengths, thermal diffusion acts to reduce thermal gradients during the 
perturbation, allowing instability even for B < Be- If the ratio of magnetic to thermal diffusivities, 
rj/lC, is small, the fluid elements retain their destabilizing magnetic buoyancy, but the stabilizing 
thermal buoyancy is washed out (e.g., Acheson 1979). The stability criterion is as given by 
equation (|36[), except that the thermal buoyancy term is reduced by a factor rj/JC. When there is 
no composition gradient, the critical magnetic field Be becomes B'^ = (r///C)^/^-Bc- We plot B'^ as 
a dotted line in Figure |5[ 

To estimate in the ocean, we use the heat equation (eq. [P) to define a thermal time 
^therm = CplP' / pK = SKCp'tj'^ / AacT^ , SO that r]/)C = tthermAdiff- The heat transport in the ocean is 
by electron conduction with K = Tr^n-e /c^Tr/Sm^,, and we use the results of ^ for tdiflf, giving 



rj _2x 10-4 11+^^2. (Z_ 



Taking Ag 



1> 514 = 1-9) Ti = 4/3 and /ig = 2, the critical magnetic field is 



B'^ = 10^ G 



1 + 



9[x, 



When x S> 1, we find 



i?^ = 5 X 10^ G y% 



.1/2 



5/8 



0.01 rhEdd/ Uo. 



(43) 



(44) 



m 



0.01 rn-Edd 



-1/2 ^^y/2 
,30/ 



(45) 



a similar scaling with depth to Be- 



To summarize, we find that instability requires B ^ lO^'' G at the base of the ocean for 
m = 0.01 rhEdd) or B ^ 10^^ G for m = 0.1 mEdd- We do not expect buoyancy instability to 
operate in the crust, as long as the rigidity of the crust is able to counteract the buoyancy force. 
The detailed outcome of the instability will depend on its non- linear development. However, it 
seems reasonable to take the critical magnetic field of equation (|4^) as an upper limit on the 
strength of any buried field at the boundary between the ocean and crust. 



7. Discussion and Conclusions 

We have investigated whether the magnetic field of a LMXB neutron star may be 
diamagnetically screened by the accreted matter. The magnetic profile in the outer layers is set 
by the competition between downwards advection and compression by accretion, and upwards 
transport by Ohmic diffusion. In steady-state, we have shown that the magnetic field strength 
decreases through the outer crust and ocean by n orders of magnitude, where n !^ m/0.02 rhEdd 
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(see §^). The steady-state profiles are stable to buoyancy instabilities provided the magnetic field 
in the ocean and crust is ^ lO^'^-lO^^ G. 

Our results are summarized in Figure ^. Arbitrarily choosing B = 10^ G at the crust /ocean 
boundary, we plot the magnetic field at the surface of the ocean for different accretion rates. 
At large accretion rates, the magnetic field strength decreases dramatically through the ocean, 
and the stellar field is screened by the accreted layer. At low accretion rates, the magnetic field 
penetrates into the accreted layer, and screening is ineffective. The dashed and dotted lines in 
Figure ^ show, with and without thermal diffusion (see §^ for the physics of this difference), the 
critical magnetic field needed at the top of the crust for magnetic buoyancy instability to occur in 
the ocean. This provides an upper limit to the strength of any buried field, since, for larger field 
strengths, we expect additional upwards transport of magnetic flux by buoyancy instabilities. 

We have shown that for accretion rates ^ 0.01 mEdd) microscopic Ohmic diffusion prevents 
screening of the neutron star magnetic field by the accreted matter. At larger accretion rates, 
field burial may be possible; however, we stress that we have considered only a plane-parallel 
model with a simple field geometry. For example, consideration of a more complex field geometry 
in a spherical shell may indicate that additional modes of magnetic instability are available. 
In addition, the burial process itself is not well- understood although some calculations of the 
spreading of matter initially confined to the polar cap of a strongly magnetized neutron star with 
B ~ lO^'^ G have been carried out (Hameury et al. 1983; Brown & Bildsten 1998; Litwin, Brown, 
& Rosner 2001). 

Our results may explain why SAX J1808. 4-3658 is the only known LMXB to show pulsations 
in its persistent emission (upper limits to the pulsed fraction in other sources are ^ 1%; Leahy, 
Eisner & Weisskopf 1983; Mereghetti & Grindlay 1987; Wood et al. 1991; Vaughan et al. 1994). 
This source has similar X-ray spectral and variability properties to other LMXBs (van der 
Klis 2000, and references therein), but has an unusually low accretion rate. The time-averaged 
accretion rate is « 10^^^ Mq yr~^ « 10^'^ "^-Edd (Chakrabarty & Morgan 1998; Bildsten & 
Chakrabarty 2001), and at the peak of outburst reaches ^ 0.03 mEdd- This accretion rate is 
similar to that at which we expect screening to become ineffective. This leads us to suggest that 
perhaps most LMXBs have their dipole magnetic fields screened by the accreted matter, whereas 
SAX J1808. 4-3658 accretes slowly enough that the magnetic field is able to rapidly penetrate the 
freshly-accreted material. The other faint Galactic center transients discovered with BeppoSAX 
(e.g.. King 2000) may therefore be promising sources to search for persistent pulsations. 

We have concentrated on the steady-state magnetic field profiles. However, there may be 
important time-dependent evolution of the magnetic field. In particular, if the magnetic field is 
screened, then we expect it to emerge once accretion ceases. The timescale for this depends on 
the depth of the screening currents (Young & Chanmugam 1995). For the inner crust, the Ohmic 
diffusion time is 10^-10^ years, always less than the accretion time (see Figure 9 of Brown &: 
Bildsten 1998). Thus we expect the relevant timescale to be the Ohmic time across the outer 
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crust, which is extremely short, only 100-1000 years. It therefore seems unlikely that screening 
alone can explain the low magnetic fields of the millisecond radio pulsars. Potentially observable 
time-dependent effects may occur in those transiently-accreting systems whose accretion rate 
during outburst is large enough for screening to be important. After an outburst, we would expect 
the magnetic field profile to relax by Ohmic diffusion on a timescale comparable to the outburst 
duration. 

Our results show that evolution of magnetic field in the outer crust is sensitive to the level of 
impurity there, as suggested by Schatz et al. (1999). For a primordial crust with a low level of 
impurity, as assumed by previous workers, phonon scattering dominates the conductivity. However, 
as shown by recent calculations of nucleosynthesis during H/He burning in the atmosphere (Schatz 
et al. 1998, 1999; Koike et al. 1999), the crust is likely very impure. In this case, the conductivity 
is set by scattering from ions, leading to shorter Ohmic dissipation times than previously assumed. 

Wc have not discussed the layer of H/He that lies above the ocean. Because of the different 
composition, the Ohmic diffusion time in this layer is longer than that in the ocean. However, we 
find that the magnetic profile in this layer is determined by a different process, thermomagnetic 
drift (Urpin & Yakovlev 1980) driven by the heat flux from nuclear burning. The evolution of the 
magnetic field in this layer is important for understanding Type I X-ray bursts. In particular, the 
drifting oscillations observed during bursts with RXTE (see Strohmaycr 1999 and van der Klis 
2000 for reviews) have been interpreted as due to a differentially rotating atmosphere. However, 
Gumming & Bildsten (2000) pointed out that a poloidal field ^ 10^ G is large enough to prevent 
differential rotation. Type I X-ray bursts are thus promising probes of the magnetic field in the 
surface layers. We address magnetic field evolution in the H/He layers in a separate paper. 
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Fig. 1. — Temperature and conductivity profiles in the ocean and outer crust. We take the 
temperature and column depth at the top of the ocean from the X-ray burst ignition models 
of Gumming &: Bildsten (2000) for rh = 0.01 (solid curves) and 0.1 rfiEdd (dotted curves), and from 
the steadily burning model of Schatz et al. (1999) for rh = rriEdd (dashed curves). We assume a 
single species of nuclei with ^ = 60 and Z = 30, and assume a constant heat flux of F = 150 keV 
per nucleon. We indicate the ocean-crust boundary with a filled circle (the mEdd model remains 
liquid for y < 10^^ g cm~^). For 0.01 and 0.1 mEdd) the two different curves beyond the filled circle 
are for either electron-ion scattering or phonon scattering in the crust. 
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Fig. 2. — Accretion time (solid line) and Ohmic diffusion time (dotted line) across a scale heiglit, 
for (top to bottom panels) m = 0.01, 0.1 and 1 mEdd- The two solutions in the crust are for phonon 
or electron-ion scattering (see Figure 1). 
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Fig. 3. — The ratio of Ohmic diffusion time across a pressure scale height tdiff to the accretion 
flow time across a pressure scale height taccD for rh = 0.01 (solid), 0.1 (dotted) and 1 rfiEdd (dashed 
line). The two solutions in the crust are for phonon or electron-ion scattering. 
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Fig. 4. — Steady-state profiles of the magnetic field, for (left to right) m = 0.006, 0.01, 0.03, 0.1, 0.3 
and 1 mEdd- We plot the magnetic field relative to the value at y = 10^^ g cm^^. The solid curves 
show the profiles for electron-ion scattering in the crust. For m = 0.01 and 0.1 mEdd) the dotted 
curves show the profiles for phonon scattering. 
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Fig. 5. — The critical magnetic field for buoyancy instability. Taking the gradient of the steady- 
state magnetic profile (shown as a dashed line, arbitrarily scaled to 10"'^'^ G at the base), we calculate 
the critical magnetic field at which buoyancy instability will occur. The solid line shows the 
critical field ignoring thermal diffusion, the dotted line includes thermal diffusion, which reduces 
the effective thermal buoyancy. For clarity, we adopt only models with electron-ion scattering in 
the crust for this figure. 
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Fig. 6. — A summary of our results. The solid lines indicate the screening of the magnetic field 
by the ocean. Assuming a value of 10^ G at the base of the ocean (upper solid line), we plot the 
B field at the ocean surface (lower solid line). The magnetic field strength decreases by roughly 
m/0.02 rn-Edd orders of magnitude through the ocean. The dashed and dotted lines show the B 
field needed at the ocean floor for magnetic buoyancy instability to occur in the ocean, with and 
without the effects of thermal diffusion, respectively. The dashed line gives an upper limit to the 
strength of any buried field. 



